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Abstract 

We establish the existence, uniqueness and main properties of the funda- 
mental solutions for the fractional porous medium equation introduced in |45j . 
They are self-similar functions of the form u(x,t) = t~ a f(\x\ t~@) with suit- 
able a and f3. As a main application of this construction, we prove that the 
asymptotic behaviour of general solutions is represented by such special solu- 
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qualitative properties of the equation we prove an Aleksandrov reflection prin- 
ciple. 
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1 Introduction 



Our main goal is to determine the existence, uniqueness and main properties of the 
solution of the equation 

(1.1) d t u + {-A) s {u m ) = 0, xeR N ,t>0, 

with fractional exponent s G (0, 1), and m > 0. This equation, or better said family of 
equations, is one of the standard models of nonlinear diffusion involving long-distance 
effects in the form of fractional Laplacian operators, which are the most representative 
nonlocal operators of elliptic type. We recall that the fractional Laplacian operator 
is a kind of differentiation operator of order 2s, for arbitrary s G (0, 1), that can be 
conveniently defined through its Fourier Transform symbol, which is |£| 2s , so that for 
s = 1 we recover the standard Laplacian. A major difference between the standard 
and the fractional Laplacian is best seen in the stochastic point of view, and it consists 
in taking into account long-range interactions in the latter, which explains features 
that we will see below, like enhanced propagation with the appearance of fat tails 
at long distances (to be compared with the typical exponentially small tails of the 
standard diffusion or the compactly supported solutions of porous medium flows). 
This is known as anomalous diffusion. There is a wide literature on the subject, both 
for its relevance to Analysis, PDEs, Potential Theory, Stochastic Processes, and for 
the growing number of applications in Mechanics and other applied fields. See basic 
information in [H EH Hfl E21 ESI [57] . 



We consider the Cauchy problem for Equation (1.1) taking as initial data a Dirac 
delta, 

(1.2) u(x,0) = M6(x) M>0. 

Solutions with such data are called fundamental solutions in the linear theory, and we 
will keep that name though their relevance is different in the nonlinear context. We 
use the concept of continuous and nonnegative weak solution introduced in [32], [46J, 
for which there is a well-developed theory when the datum is a function in L 1 (M. N ). 

Putting s = 1 we recover the standard Porous Medium Equation (PME), where the 
question under discussion is well-known, see the comments in the next section. The 
exponent m varies in principle in the range m > 1, but the methods extend to the 
linear case m — 1. and even to the fast diffusion range (FDE) m < 1 on the condition 
that m > m c = max{(iV — 2s) /N, 0}. Such type of restriction on m from below 
carries over from the PME-FDE theory, 
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A main result of this paper consists in proving that there exists a unique fundamental 
solution of Problem (!.!)-( 1.2) as follows. 



Theorem 1.1 For every choice of parameters s 6 (0, 1) and m > m c = max{(iV — 



2s)/N, 0}, and every M > 0, equation (1.1) admits a unique fundamental solution; 



it is a nonnegative and continuous weak solution for t > and takes the initial data 



(1.2) as a trace in the sense of Radon measures. Such solution has the self-similar 



form 

(1.3) u* M (x,t) = r a f(\x\r?) 

for suitable a and (3 that can be calculated in terms of N and s in a dimensional way, 
precisely 

N 1 

!- 4 a = Tu 7^^> P 



N{m-l) + 2s , N(m-l)+2s 

The profile function f(r), r > 0, is a bounded and Holder continuous function^ it is 
positive everywhere, it is monotone and goes to zero at infinity. 

More precisely, the initial data are taken in the weak sense of measures 

(1.5) lim / u(x,t)4>(x)dx = M<f>{0) 

Jr n 

for all cf) G C&(R ), the space of continuous and bounded functions in M. N . We will 



call these self-similar solutions of Problem ( 1 .1 )-( 1.2 ) with given M > the Baren- 
blatt solutions of the fractional diffusion model by analogy with the PME and other 
prominent studied cases. The form of the exponents explains the already mentioned 
restriction on m from below. The result is proved in Sections [3] to [9| 

We then prove that the asymptotic behaviour as t — > oo of the class of nonnegative 



weak solutions of (1.1) with finite mass (i.e., Ju(x,t)dx < oo) is given in first 



approximation by the family u* M {x } t), i. e., the Barenblatt solutions are the attractors 



in that class of solutions. See Section flOl Theorem 10.1 



In the limit M — > oo we obtain a special solution with a fixed isolated singularity 
at x = that has separate- variables form and we call very singular solution (VSS) 
by analogy with the standard Fast Diffusion Equation. This happens precisely in the 



/is actually a smooth function, but such result is not discussed or used here. 
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range m c < m < mi = N/ (N + 2s), see Theorems 13.1 , 13.2 The result represents a 



marked difference with the case s — 1 where VSS exist in the larger range m c < m < 1. 
It is another manifestation of the long-range interactions of the fractional Laplacian, 
that avoids some of the purely local estimates of the standard FDE with the classical 
Laplacian operator; indeed, extrapolation of the standard estimates would justify the 
existence of a VSS in cases there is none. Let us recall that the class of VSS has 
played a role in the study of nonlinear parabolic equations since the seminal paper of 
Brezis, Peletier and Terman [201. 



We then devote one Section to constructing eternal solutions in the critical exponent 
case m = (N — 2s) /N. They have exponential self-similarity. 

The study of evolution equations with fractional Laplacian operators implies the 
need to develop in this setting tools that have been successful in the standard and 
parabolic theory. One of them is the Aleksandrov reflection principle that we state 
and prove in the elliptic and parabolic settings in Section [15} 

A number of appendices and a section of comments and extensions close the paper. 



2 Motivation and historical perspective 

The question of finding fundamental solutions of elliptic and parabolic equations is 
one of paramount importance in the study of linear elliptic and parabolic equations, 
see e. g. [28> 32]. In the nonlinear theory their role is not so apparent, but it has been 
shown that they can be an important tool in the existence and regularity theory, and 
very important in the description of the asymptotic behaviour. 

Generalizing the classical heat equation, a prominent case at hand in the evolution 
theories is the Porous Medium Equation (PME), d t u — A(u m ) — 0, m > 1, introduced 
in the last century in connection with a number of physical applications and exten- 
sively studied as a prototype of nonlinear diffusive evolution with interesting analysis 
and geometry (free boundaries). Fundamental solutions were discovered in the 1950's 
by Zeldovich and Kompanyeets [59] and Barenblatt [5], and later by Pattle [47J, the 
latter made a complete description. They have been usually called Barenblatt solu- 
tions in the literature. This discovery was so to say the starting point of the rigorous 
mathematical theory that has been gradually developed since then. Their role in the 
asymptotic behaviour of general solutions of the PME was established in papers by 
Kamin [351 136] , Friedman and Kamin [29] and the author [51] . They have also played 
an important role in the existence and regularity theory in the work by many authors. 
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The monograph [56] contains a detailed exposition of the topic. 

The surprising relation between existence and uniqueness of fundamental solutions 
and precise asymptotic behaviour relies on the existence of a scaling group under 
which the solutions of the PME are invariant. It implies that a fundamental solution 
is in fact a self-similar function. This is what we call a Barenblatt solution. We will 
see below how the scaling group works in the case of equation (1.1). Self-similarity 
plays a big role in our understanding of fundamental processes in mathematics and 
mechanics, as explained by Prof. Barenblatt in his books [8] [9] . 

Following the analysis of Barenblatt solutions for the PME, other equations have been 
explored. To name a few, the Fast Diffusion Equation, which is dtu — A(u m ) = 0, 
with m < 1. Though it looks formally the same, the qualitative theory changes. The 
role of the Barenblatt solutions remains basically unchanged only for m > (N — 2)/N, 
while for m < (N — 2)/N the whole functional setting changes abruptly, cf. [55J and 
[T5J [T7|. Barenblatt solutions have been constructed for the p-Laplacian equation, 
d t u — A p (u) = |6]; in this case they exist in the range 2N/(N + 1) < p < oo and 
their role and properties are quite similar (loosely speaking) to the porous medium 
case. These solutions have been used extensively in the regularity theory [26] and in 
the asymptotic behaviour [2D]. The ideas can be extended to the doubly- nonlinear 
diffusion equation d t u — A p {u m ) = 0, as proposed in [BJ, see [2] for recent results, 
and to other models, not necessarily of diffusive type. Among such models we count 
conservation laws, [44J. 

In some cases the fundamental solutions have similarity exponents that cannot be 
calculated from dimensional considerations, giving rise to so-called anomalous ex- 
ponents, it is also described as self-similarity of the second kind. Then the initial 
distribution is not a Dirac delta but a more general mass distribution, in other words 
some singularity located at x — 0, t — 0. We say that it is a solution with an isolated 
singularity, or a source-type solution, but not a fundamental solution. There is a 
wide literature on solutions with isolated singularities, both in elliptic and parabolic 
problems, among them particular attention has given to the so-called very singular 
solutions, cf . [TTJ [2DJ |37J |38j [58] • Again, the influence of Prof. Barenblatt has been 
felt in the form of interesting models from the physical sciences. One of them is the 
Barenblatt equation of elasto-plastic filtration, proposed in [5] , where the source-type 
solutions have been studied in [TDJ|39], see also [33J. They correspond to self-similarity 
of the second kind, with anomalous exponents. Another model is the equation of tur- 
bulent bursts [7] treated mathematically in [UJ, [30J and [31], where self-similarity 
with anomalous exponents is also found as asymptotic behaviour of a wide class of 
solutions. Anomalous exponents are discussed in [3J in connection with focussing 



6 



problems, and they are found in the Fast Diffusion Equation under the critical expo- 
nent [12, HH [55], but in these cases they do not correspond to source-type solutions. 

Nonlinear diffusion with nonlocal operators. There are many other in- 
stances of the paradigm we are discussing. Recently there has been a remarkable 
interest in mathematical models of diffusion involving long-range effects represented 
by fractional Laplacian operators and other singular integral operators. This responds 
to a serious motivation from the physical sciences and has produced so far quite in- 
teresting mathematical developments. Let us add some comments to the information 
that opened the paper on fractional Laplacian operators. For functional analysis con- 
siderations, it is usually best to consider the inverse operator (— A) _s , which happens 
to be the integral operator associated to the Riesz kernels, as we will do often in this 
paper. Also, a special interest has the case N = 1 where the operator (— A) s is the 
natural realization of the concept of (d x ) 2s (up to a constant factor) in the form of a 
positive symmetric operator. 

Regarding nonlinear evolution models, the author has been involved in the analysis 
of two of such models, called fractional porous medium equations. One of them is 



Equation (1.1) that we study here, while the other model uses the equation 



(2.1) d t u - V • (uV(-A)-V) = 0, < s < 1 . 

Both models have quite different properties, see a survey of recent results in [57]. 
The construction of the Barenblatt solutions for this latter model has been performed 
by Caffarelli and the author in [22]. The asymptotic behaviour is obtained in [22] 
using sophisticated entropy and obstacle problem methods. An explicit form for the 
Barenblatt solutions is found in 



(2.2) U(x,t) = Cir Q (l - c 2 \x\H- 2a/N ) 1 - s with a = N/{N + 2 - 2s). 

In the limit s — > 1 we obtain quite interesting fundamental solutions of the non- 
diffusive limit Ut — V • («V(-A) _1 m) = 0, of interest in superconductivity and super- 
fluidity. The solutions take the form 

(2.3) u* M (x,t) = j X B R (t)(x), R(t) = ct 1 ' N , c>0, 

as described in [13] and [50]. Each one represents a round vortex patch that expands 
to fill the space as t — > oo. 



The linear case. Solving Equation (1.1) in the linear case m — 1 is easier. As 
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explained in [36], when m — 1 and < s < 1 the evolution has the integral represen- 
tation 



(2.4) u(x,t) = / K s (x — z,t)u (z) dz, 



where K s has Fourier transform K s (^,t) = e~'^' 2s *. This means that, for < s < 1, 
the kernel K s has the form K s (x,t) = t~ N / 2s F(\x\ t~ l l 2s ) for some profile function F 
that is positive and decreasing, and behaves at infinity like F(r) ~ r -{N+2s)^ gee ^gj 
When s = 1/2, the kernel is explicit, 

(2.5) K l/2 (x,t) = C N t(x 2 + t 2 p N+i y 2 . 

If s = 1 the function Ki(x, t) is the Gaussian heat kernel. The linear model has been 
well studied by probabilists. 



However, an integral representation of the evolution like (2.4) is not available in 
the nonlinear case; the tools to treat the nonlinear case were developed in [35], [3E]. 
Below, we perform the analysis of existence, uniqueness, properties and applications 



of Barenblatt solutions for equation (1.1). We will concentrate on the difficulties 



caused by the nonlinear form of the equation. Our Barenblatt solutions will play the 
role of the kernel K s in the study of asymptotic behaviour. Our results below include 
the linear case as a particular instance. Let us recall that the nonlinear techniques 
are quite different from the linear ones. 



3 Preliminaries 

Before we start the construction we need to review some basic facts. 



Scaling. Equation (1.1) is invariant under translations in space and time, and also 
under the scaling group, and this fact will play a big role in what follows. Let first 
observe that for every solution u(x, t) and constants A,B,C > the function 

(3.1) u(x,t) = Au(Bx,Ct) 

is again a solution of (1.1) if C — A m ~ l B 2s . This generates the whole scaling group. 



If moreover, we impose the condition that the solutions have constant finite mass 
(i. e., integral in space), then A = B N , which implies that the group is reduced to the 
one-parameter family u(x, t) = B N u(Bx, B N(jri ~ l ^ +2s t), which is usually written in 
terms of A = B N ^ m ^ +2a as 

(3.2) (T x u) (x, t) = \ a u(\ p x, X t) , 
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with scaling exponents given by formula (1.4): 



a = — — — — , /3 



N{m - 1) + 2s N(m - 1) + 2s ' 

Note that a and /3 > iff m > m c = (N — 2s) /N. The values of both parameters 
will be fixed in what follows. 

Another application of scaling consists in reducing solutions to unit mass. Indeed, if 
u(x, t) is a solution with J u(x, 0) dx = M > then 

(3.3) u(x, t) = M~ l u{x, M- {m - l) t), 
is another solution with unit initial mass, J u(x, 0) dx = 1. 

Potential equation. Take the convolution U(x,t) = u(x,t) *Z 2s (x), where X 2s is 
the Riesz kernel 

(3.4) l 2s (x) = C N , s \x\~ (N ~ 2s) , 

that is the kernel of operator (— A)~ s , < s < lj^] Use of this kernel restricts the 
range of application toO < s < 1/2 when N — 1. This is a typical difficulty of 
potential theory that is found in the standard heat and porous medium equations in 
dimensions = 1,2. We will discuss the case N = 1, 1/2 < s < 1, separately in 
Section |9l 

After this caveat, we resume the theory. We apply the kernel to a solution u(x,t): 

Then (— A) S U = u, and using the equation we get equation (PE): 

(3.5) U t = ((-A)- s [/) t = -(-A)- S (-A) s u m = -u 



in 



This is easy to justify since u m (t) G L 2 with u t = (—A) s u m G L 1 for every t > 0, 
according to the theory. 

The scaling group acts also on the solutions of the PE: If A, B > are positive and 
U is a solution of the PE, then so is 



(3.6) U(x, t) = A U(Bx, Ct), with C = A m - l B 2m 



2 For the value of Cn, s see formula ( |15.2[ ) below. 



is again a solution of the same equation with suitably rescaled initial data. 

The idea of using the potential equation to prove uniqueness of solutions of diffusive 
equations with measure data goes back to Pierre's work for the standard PME |49j. 
This has been followed later by a number of authors in various contexts. 



4 Existence of solutions with measure data 

Here is the basic result in this topic. «M + (R JV ) is the space of bounded and nonneg- 
ative Radon measures on M. N . 

Theorem 4.1 For every fi e A^ + (M Ar ) there exists a nonnegative and continuous 
weak solution of Equation (1.1) in Q = IR x (0, oo) taking initial data fi in the sense 



that for every ip G C|(R ) we have 

(4.1) lim / u(x,t)(p(x) dx = I (p(x)dp(x) 



Proof, (i) Existence comes from approximation of the initial data \i with a smooth 
mollifier sequence p e [x\ thus we take smooth initial data fi e = fi * p £ . Using the 
results of jl6] we know that the corresponding solutions u £ exist and are uniformly 
bounded in L°°(0,oo : L 1 (R N )) and L°°(R N x (r, oo)). The decay rate with time, 
||w(t)||oo < Ct~ a , is a main tool, and it holds uniformly for the whole sequence. They 
also have uniform bounds on the energy J f °° f |(— A) s ^ 2 u m \ 2 dxdt, cf. [36], formula 
(8.5). 

(ii) We pass to the limit e — > 0. Convergence holds in principle up to subsequences 
e n — > 0. The solutions u £ and the limit u are uniformly C a in space and time for all 
t > t > 0. The limit is easily proved to be a weak solution of the equation. 

(iii) The initial data are taken in the sense of initial traces, since the weak formulation 
passes to the limit for t > r > and the second integral term is uniformly small for 
small t e (0, r): if L = (-A) s then 

u m L V dxdt\ < I |Ks)||™ _1 ( f u\L(p\dx] ds. 



Now we use that \L(p\ < C, Judx < C and the decay ^(t)^ < Ct a with a 
N/{N{m - 1) + 2s). Therefore, 

4 -2 aim- 1 = — v - ' < 1 

v ' K ' N(m-l) + 2s 
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and the integral in time converges at t — 0. This calculation is valid for the u e and 
also for the limit u. See a similar argument in [50J, Section 3. So we have initial 
traces in the sense of distributions. Even better, since the test functions can be taken 
in W 2s '°°(R N ) for instance. 

(iv) The previous argument as t — > works for m > 1. For m c < m < 1, we use good 
test functions and Holder estimates: 

(4.3) | jj Lu m <fdxdt\ < (^jj udxdt^j (^jj (L<p) 1/{1 - m ^ <Ct x - m . 

□ 



5 Approximate data for the fundamental solution 

1. In our study of the solutions with Dirac deltas as initial data, we will use some 
monotone approximations. The simplest approach is to take a smooth convolution 
kernel p G C^°{R N ), p > 0, J pdx = 1. We also take p to be radial, i.e., a function 
of \x\ and supported in the ball of radius 1. Put p e = e~ N p(x/e), and define the 
regularization V e = p e * V of the Riesz kernel V(x) = l2 S ( 



where C\ is the constant C(N, s) of (3.4). Note the scaling 



V £ ,(x) = V £ (x/p)/p N ~ 2s 

that allows to reduce the calculations to the case e = 1. Clearly V e — > V in L; 1 0C (IR 7V ). 
We want to examine further the approximation of V by V e . 

• We claim that 

V e {x) > V(x) V|x| > e. 

Proof: It is enough to prove it for e — 1. Using spherical coordinates we write y = ra, 
\a = 1|, Jj a ,| <1 dx = un, and then 

V 1 (x) = C 1 [ p{r)dr [ l - w ^r r n - l da = Nu N d f p{r)r N - l ${r)dr 



J\y\=r \ x il 
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where $(r) := jv y \ =r \x — y\ N+2s da. So we are reduced to study the properties of 
last integral. This is where we use the fact that for z ^ 

AV(z) = (N- 2s) (2 - 2s)V(z)/\z\ 2 > 0, 

so the integrand function is subharmonic in B r (x) if r < \x\. We conclude that 

J\y\=r \ x y\ 

and the result is true. In fact, Evans [28], page 26, proves that 

$'( r ) = (r/ N ) / AV(x -y)dy=-^ 1 [ AV(x - y) dy 

JB r (0) r J B r {0) 

hence $'(r) > for all r < \x\. 

• We now estimate the error for large |x| > 2. In that situation we have 
$'(r) < c 2 r\AV(x)\, $(r) - $(0) < c 3 r 2 V(x)\x\~ 2 , 

so that 

i \ 



V^x) - V(x) < c 4 ( J p(r)r N+1 drj V(x)\x\~ 2 = c 5 V(x)\x 

Using now the rescaling rule, we get interesting results for V e . 

Proposition 5.1 There is a constant K > depending on N and s such that 

e 2 



(5.1) V(x) < V £ (x) < V(x) [1 + K 
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\X\ 

holds for all e > and all \x\ > 2e. There is another constant C such that 

(5.2) V £ (x) < CV(x) 

for alle>0 and all x. Finally, V e - V -»• in L\R N ). 
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For the last part, note that for x away from zero |T4(x) — V(x)| < Ce 2 \x\~ N ~ 2+2s 
which is integrable at infinity, hence the L 1 norm in that region is small if e — > 0. 
Near zero \V £ (x) — V(x)\ < CV(x) which is also integrable, hence uniformly small in 
a small ball. 

2. The following construction supplies an alternative approach that may also be 
useful. We consider V n (x) = inf{V(x),n} for n > 1. Note that this is just a rescaling 
on Vi(x) of the form V n (x) = n V^n 1 /^ -2 ^). Clearly, V n (x) increases monotonically 
to V(x) as n ->■ oo, V n - V ->■ in L l {R. N ). 

Proposition 5.2 Consider the function <f) n (x) := (— A) S V^(^)- Then, j <f> n (x) dx = 
§(—A) s V dx = 1 and <f> n (x) > everywhere in M. N . The sequence <p n = (— A) s V n (x) 
is a positive approximation sequence to the Dirac delta. 

Proof, (i) We concentrate on <f> 1 . On one hand, V\ has a flat plateau in a ball \x\ < R 
where it attains the maximum value 1, and by a simple calculation with the explicit 
formula for the operator we get (— A) s Vi(x) > 0, and <p\ is bounded and continuous 
for \x\ < R. In order to calculate the value of (— A) s Vi(x) at the points \x\ > R where 
V±(x) = V(x) we define 

W{x) = V(x) - Vi(x) = (V(x) - 1)+. 

It is also easy to see that (— A) S W < at the points \x\ > R since W attains there 
its minimum, W(x) = 0. The integral defining (— A) S W is bounded for \x\ > R and 
there is an asymptotic estimate |(— A) S 1T^| ~ C|a;| _ ( n+2s ) as \x\ — > oo. Finally, by 
using integration with respect to a suitable cutoff function we prove that (—A) S W is 
integrable and J (— A) s W(x) dx = 0. 

Hence, (f>i(x) = (— A) s Vi(x) > also at those points. We also get the estimate 

Mx) ~ C{\x\- {n+2s) ) as \x\ -)■ oo. 

Besides, 0i(|x|) is continuous, bounded and smooth but for a possible asymptote at 
\x\ = R. Also, / (j)i{x)dx = f(-A) s Vdx = I. 

(ii) Consider now (f> n (x) := (— A) s V n (x). It follows that <p n is positive everywhere in 
R N , that n (x) = n^^-^faixn 1 /^-^). This means that / <f> n (x) dx = f<f> 1 (x)dx 
and thus <p n is a suitable approximation of the Dirac delta. □ 



13 



6 Special construction of a fundamental solution 



(i) In the existence part we use approximation of the special initial data M8(x) with 
a smooth mollifier sequence of the form Mp e (x) = Me~ N p(x/e). Without loss of 
regularity we fix M = 1 here (using mass scaling, formula ( |3.3 )). It is then clear 



that the scaling group, formula (3.2), transforms the solution u £ with data p £ into the 



solution T\u £ with data p\ £ . By uniqueness of solutions we conclude that 

(6.1) T x u £ = u Xe • 
This will be used below. 

(ii) The Potential Equation, Continuity. The corresponding solutions U £ of 
the PE (with initial data V £ ) are regular in space, since (— A) S U £ = u £ G C". Since 

(6.2) d t U £ = -u™ 

they are also C 1,a in time uniformly for t > r > 0. They are also monotone in time 
for all t > 0. 

Before examining the limit, we will establish another important property, the con- 
tinuity of the evolution orbit U(-,t) in L 1 (IR iV ). In fact, for every < t' < t we 
have 



/ 



\U{t') - U(t)\dx < j J \U t \ dxdt = j J \u m \ dxdt < oo. 



The worst case happens for t' = and we have proved above that the last integral 
can be estimated as 0(t 2sa ), that goes to zero as t — > 0. This calculation is valid for 
the approximate solutions U £ , that are classical, and also for their limits. 



We have assumed m > 1. For m < 1 using an argument like (4.3) we have 
(6.3) f(U(x,0) - U(x,t))(pdx < Ct. 

So in this case we have a uniform continuity control in L} oc (K N ). The rest is similar. 

(iii) Convergence for the Potential Equation. We have proved the lemma 
that says that U e (x, 0) < CV(x), therefore, CV(x) is a uniform upper bound for the 



initial data, and in view of (6.2) also for the solutions. 

In view of the internal regularity, we can pass to the limit along a subsequence e n — > 
to obtain a limit U*(x,t) that is a solution of the potential equation, < U* < CV. 
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As for the initial trace, in view of the convergence V £ (x) = U e (x, 0) to V in L x ( 



(Proposition 5.1 ), and the uniform continuity of the orbits U e (t) in L 1 (IR 7V ), the limit 
has also continuous time-increments in L 1 , and takes the initial data V in the sense 
that 

(6.4) \\U*(t) -y||x ^0 as t 0. 

Since d t U* < 0, U*(x,t) converges monotonically to V(x) as t — > for fixed x. 

It is clear that for positive times the following holds: (— A) s U*(x, t) must coincide 
with u, the limit of u E (x,t) = (—A) s U £ (x,t) along e n (in other words, the operator is 
closed). Since u(t) — > 5 in V this also implies that U*(t) — > V in V, i.e., U*(x, 0) = 
V(x). 

The fact that the limit U*(x,t) takes the initial data with uniform convergence 
away from zero is a consequence of the monotone convergence theorem for continuous 
functions. Around zero the convergence holds in the Marcinkiewicz space M P (M. N ), 
p = N/(N - 2s), hence in L q (R N ), 1 < q < p. 

Besides, the convergence as e — > is uniform for t > r > 0. 

We have concluded that U* is a solution of the PE corresponding to a fundamental 



solution of the original equation (1.1), and it takes the initial data V(x) in a strong 



way. We now have to solve the question of uniqueness of this limit U*. 



7 A uniqueness result for fundamental solutions 

Theorem 7.1 The fundamental solution U(x, t) of the potential equation with initial 
data V(x) is unique under the assumptions that U > 0, U is continuous in Q = 
M> N x (0, oo), and has only an initial singularity at (0,0) in the sense that U is 
bounded for t > r > 0, is continuous as t — » zero for all x ^ 0, and U(x,t) — > oo 
as (x,t) -> (0,0). 

Proof. Take two solutions U\ and Ui of the potential equation with the same initial 
data V(x). Consider U3 = U\ + 5 and consider the approximate solution Il^x^t) = 
U2(x,t + r) for a small r > 0. Then U3 is strictly larger than the continuous and 
bounded function U2 at t — 0. The first point of contact cannot happen near infinity 
because of the uniform decay of U\ and U 2 that implies strictly separation for U 3 and 
U 2 . 
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At any contact point (xq, to), to > 0, and since C/3 — U 2 > up to that time, we have 

(-A) S (U 3 - U 2 ) < 0, i.e. Ul <u 2 , 

using the representation formula for (— A) s . On the other hand, at this point 

> dt(U 3 - U2) = -u? + u™, < > v%. 

This allows to arrive at a contradiction. We conclude that Ui(x, t) + 5 > U 2 (x,t + r) 
for all 5, r > 0, hence in the limit 

Ui(x,t) >U 2 (x,t) for all (x,t) e Q . 

This concludes the proof, after reversing the roles of U\ and U 2 - 

We now know that we do not need to take subsequences in the limit as t — > 
of Section |6l Note that we could have restricted the existence time and consider 
Qt = ^ N x (0, T) instead of Q, and then the conclusion is valid for < t < T. 



8 Self-similarity of the Barenblatt solution 

We conclude here the construction of the Barenblatt solution and the proof of Theo- 
rem [Tj] for N > 1. 

• Let us settle first the question of self-similarity. Indeed, the initial function of 
the solution of the PE with data Uq(x) = C\x\~^ N ~ 2s ^ is invariant under the scaling 
T defined by (TUo){x) = AUq(Bx) if we choose A = B N ~ 2s . If moreover we have 



uniqueness of such solution, then the scaling group introduced in (3.6) will imply that 



TU(x, t) = B N - 2s U(Bx, At) with A = A m ~ 1 B 2sm = B N{m - 1)+2s 

will be a solution with the same initial data. By the claimed uniqueness TU = U, 
hence for all B > 0, x 6 M. N and t > we have 

U(x,t) = \ ( - N - 2s ^U(\Px,\t), 

where we have used the value for /3 given in the Introduction. Fix now t\ > and let 
A = l/t x . We get 

U(x,t)=t[ N - 2s)f} U(xtf J l). 
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Of course, t\ is arbitrary and can be replaced by t. Calling now U(x, 1) = H(x) we 
get 

U{x,t) =t^ N ~ 2s ^ H{xr p ). 
Applying operator (— A) s we get self-similarity for the corresponding fundamental 



solution of ( 1.1 ): 



u*(x,t) = t~ {N - 2s ^ (-A) s x (H(xr^)) = t- a F(xT p ) . 

This is a rather classical argument in the study of self-similarity. 

• End of proof of the Theorem. The proof that u*(x,t) takes the data in the sense 
of measures (and not only as distributions) follows easy from self-similarity. The 
fact that F(r) is monotone nonincreasing in r is a consequence of the Aleksandrov 



reflection principle proved in Section 15 We know by qualitative theory that F is 



always positive since u is. Same argument for the Holder continuity, n 

8.1 Properties of the profile F 

Equation. The self-similar profile F satisfies an elliptic equation 

(8.1) (— A) s F m = aF + (3y ■ VF = /3V • (yF) 
so that, putting s' = 1 — s we have 

V{-A)~ s 'F m = -(3yF, 
which in radial coordinates gives 

/■oo 

L s ,F m (r) = /3 rF(r)dr, 

J r 

where L s > the radial expression of operator (— A)~ s '. 

Dependence on the mass. The scaling group acts on the profiles Fm{t) for different 
masses M > and indeed we have 

(8.2) F M (r) =/i 2s F 1 (/i 1 - m r), M = ^ l N ( m - 1 )+ 2s . 

which reduces all calculations to the case M = 1. Since N(m — 1) + 2s > for 
M > m c we get Fm(0) — > oo as M — > oo. For m > 1 the same result happens for 
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Figure 1: Computed Barenblatt profiles for m — 1,2, 10 with s = 1/2 
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all r > 0, limM-s>oo Fm{t) = oo. However, this last limit may be finite for m < 1, see 
Section [13] to understand when and why. 



Decay at infinity. First estimate. The precise behaviour of the fundamental 
profiles F(y) = F mt N yS {y) as y — > oo is a very important question in the qualitative 
theory. It is known in the linear case m — 1, since F is given by a linear kernel K 
that decays like \y\~^ N+2s \ [IS]- The exact rate of decay for m ^ 1 is a nontrivial 
issue that we discuss next. 

To begin with, the fact that F is monotone as a function of r and also integrable in 
M N , implies that there is a constant C = C(F) > such that 

(8.3) F(r)<Cr~ N . 

We will use the same letter C for different positive constants as long as their value is 



not important in the context. We will continue with this issue in Section 12 



9 Peculiarities of one-dimensional flow 

|9[l. We add here the comments that are needed to close the case N = 1, in the 
range 1/2 < s < 1. In that case the kernel is a function that grows at infinity and 
the potential approach cannot be used in the direct way we have done before. There 
are at least 3 natural ways of addressing the difficulty 

One of them is to consider in a first stage the problem posed in a bounded domain 
Q with nice boundary, say a ball Br. Then the kernel is replaced by the Green 
function (with zero outside conditions) and this avoids considering the divergence of 
the kernel as \x\ — > oo. Once the fundamental solution for this problem is constructed, 
one passes to the limit as R — > oo and then proves uniqueness. 

A second approach is to replace the fractional operator (— A) s by a coercive operator 
like L £ = (— A)~ s + el, whose inverse has a nice kernel (the Bessel kernel Gz 8 (x)). If 
U solves the potential equation as before and L e U = u, then u solves u t = L £ U t = 
-L £ (u m ),i.e., 

(9.1) u t + (-A)V™ + eu m = 0. 

The plan is to construct a unique fundamental solution for this equation, and then 
pass to the limit and prove uniqueness. 
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The third option is to integrate in x, v(x, t) = f_ u(y, t) dy. Then v solves 
(9.2) «t = %D, L = d x (-A) s - 1 . 

THis is a kind of fractional p-Laplacian operator. 

[9[2. Let us develop the second method. We take e > and solve the equation 
u t + (-A) s u m + eu m = for x G R and t > with initial data u G L 1 (R) by just 
copying the method used in paper [46j when e = (since the new term is dissipative, 
the needed estimates still hold). We obtain a solution with similar properties (except 
for the mass conservation). Again, for fixed e > the maps S £ (t) : uq u £ (t) are 
ordered contractions in L X (IR). If uq > so is the solution u £ for all times and the 
family of solutions {u £ : e > 0} is monotone increasing as e goes down to zero. In 
this way the standard solution u of the problem u t + (— A) s (u m ) = 0, u(x, 0) = u is 
obtained in the limit e — > 0. 

We now introduce the potential functions U e using the Bessel potential instead of 
the Riesz potential, so that U e (-,t) = ((— A) s + el)~ l u £ {-, t). The U £ will satisfy, all 
of them, the potential equation U t = —u m . We prove the existence and uniqueness 
of a solution with data G2 S {x) = ((— A) s + eI)~ 1 M5(x) much as above. This implies 
that the fundamental solution UM, £ ( x ,t) of equation u t + (—A) s u m + eu m = exists 
and is unique. We leave these details to the reader. 

It follows easily that Mjv/, £ (x,t) is monotone increasing as e goes down to zero and 



the limit is a fundamental solution of Equation (1.1) for N — 1. We then prove that 
the limit is a minimal element among such possible fundamental solutions, and then 
conservation of mass implies that the fundamental solution must be unique. The 
result of Section [8] follows then by easy changes in the argument. 

Note. We have collected some properties of the Bessel potentials in an Appendix 
at the end of the paper for the reader's convenience. Note that these potentials have 
exponential decay as |x| — > oo. Also (in N = 1) they are singular unbounded at x = 
for 2s = a < 1, but they are bounded for 2 s > 1. 



10 Asymptotic behaviour of general solutions 

We use the methods of the monograph [52] to get the following general theorem. We 
recall that all our solutions are nonnegative. 
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Theorem 10.1 Let u = n G M+(R N ), let M = fi(R N ) and let u* M be the self- 
similar Barenblatt solution with mass M . Then we have 

(10.1) lim t a \u(x,t) -u* M (x,t;M)\ = 0, 

t— >oo 

and the convergence is uniform in R N . 

Proof. I. Firstly, we make the proof under the assumption that d/i = u (x) dx where 
Uq is a bounded function with compact support. We divide this proof into several 
steps. 

(i) We may assume that there are K > big enough and r > such that u (x) < 
u* k (x,t). The Maximum Principle implies then that 

u(x, t) < u* K (x, t + t) for all xeR N , t> 0.. 

We now perform the mass-preserving scaling transformation 

(T x u)(x,t) = \ a u(x\P,\t). 

It is easy to see that we have a similar estimate for the family u\ = T\u: 

(T x u)(x,t) < u* K (x,t + (r/A)) forall xeR N ,t>0. 

The family (T x u) is bounded in Z^nLTl C a in any set of the form Q = B R (0) x (1/2, 2) 
hence it converges to a solution u of the equation for t > 0. The limit satisfies the 
bound u(x,t) < u* K (x,t) for all x G R N , t > 0. 

(ii) We claim that u takes the initial data MS in the sense of measures. Indeed, away 
from zero the convergence of the limit u(x, t) to zero as t — > is uniform, because it is 
bounded above by a tail of the self-similar Barenblatt solution u* K (x, t). This leaves 
as ppssible initial trace a Dirac delta. It only remains to recall that the total mass of 
all the family (T\u) is the same, M, at all times. 

(iii) At this point we pass to the corresponding solutions U\ = (A)~ s (T A w). They 
are solutions of the PE. The limit U is also a solution of the PE with the already 
mentioned regularity and monotonicity properties. The initial datum of U (x, t) is M 
times V(x) in a weak sense, since 

f(V W - U( X , WW * ^ M(-Ar>(0) -/«*,«)<-*)-«*) * , 
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which goes to zero since ((x) = (—A) s <fi(x) is an acceptable test function for the 
convergence in the sense of measures of u(-,t) to M 5. 

On the other hand, the family U(t) converges in L 1 (]R 7V ) and monotonically to the 
initial trace that is M V(x), as we have shown. The properties of the uniqueness the- 



orem (7.1) are met, so that we conclude that U is the fundamental solution U^ix^t). 
(iv) Applying (— A) s to U = U* we obtain the consequence u = u* M . We thus have 

(T\u)(x, 1) —7- u* M (x, 1) uniformly in x G M. N . 



It is just a routine calculation to transform this expression into the convergence (10.1 ), 
see for instance [SB], chapter 18. 

II. For a general initial function w L 1 (M Ar ) we have to work a bit more. First we 
fix an 5 > and truncate uq above and near infinity to fall into the previous case, 
namely 

< u QjS (x) < u (x), \\uo,s\U = M 5 > M - 5. 

In this way we have uniform convergence of the rescaled versions Ug(x, 1) to u* M (x, 1; M$ 
at t — 1. Since the sequence {(T\Us)(x, 1) : 5 > 0} is uniformly bounded and mono- 
tone in 5 we get the convergence of the {(T\Us)(x, 1) : 5 > 0} to u* M (x,t; M) in L 1 
and since they are bounded, they are uniformly C a , hence convergence is uniform on 
any fixed ball. 

In order to examine the convergence of (T\Us)(x, 1) outside a big ball we first observe 
that all the functions are bounded in L 1 D L°°, hence they are C a with a uniform 
constant and exponent, which makes them a compact set of functions, hence we can 
extract a uniformly convergent subsequence. On the other hand, the L 1 norm of of 
(T\u)(x, 1) is small in the outer domain as a consequence of the following clever trick: 
fixing 5 > small, such L 1 norm is bounded by 5 plus the norm of of (T\iis)(x, 1) , 
and this one is controlled in that domain by the tail of u* K (x, 1 + r/A), with K = K{8) 
large. By interpolation between the L 1 and C ia norms we get a small norm in L 
uniformly in 5 for e small. 



oo 

) 



(iii) For a measure as initial data, we displace the origin of times to t = r > and 
may then assume that uq is integrable and bounded. □ 

• We have a strong version of uniqueness of self-similar solution, using the asymptotic 
behaviour: 



Corollary 10.2 Let u(x,t) be a nonnegative weak solution of Equation (1.1), and 



assume that u is self-similar and the mass if finite, i.e., J u(x,t)dx < oo for t > 
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0. Then, u is one of the Barenblatt solutions mentioned in the previous theorem, 
u(x,t) = UM(x,t), where M = J u(x,t) dx, which is constant in time. 



11 Nonexistence for small m < m c 

When we go below the critical exponent m c , fundamental solutions cease to exist. This 
was proved for the standard Fast Diffusion Equation (s = 1, < m < (N — 2)/N) 
by Brezis and Friedman [Tj5]. Without entering in full details on the issue, we give 
here a simple proof of nonexistence in the fractional case, based on the existence 



of the scaling group (3.2). The main point to take into account is that when m < 



m c the similarity exponent /3 of formula (1.4) is negative, i.e., the transformation 



involves space contraction instead of expansion. Therefore, if we consider a sequence 
of integrable functions approximating the Dirac delta of the form 

Uon(x) = n N Uoi(nx), with J Uoi(x) dx = 1, 

and say u 0n G C^°{R N ), then the sequence of solutions satisfy 

u n (x, t) = n N Ui(nx, n l ^t) . 

Passing to the limit n — > oo for fixed t > the formula implies that u n (-,t) — > 5 
(thanks to the fact that 0), and we conclude two things: first, that we do not 

find in the limit the expected fundamental solution, and second, that the result can 
be interpreted as saying that, under this evolution equation, an initial Dirac delta 
does not spread with time and the "physical solution" is just u* M (x,t) = M5(x). 

In the critical case m = m c we have 1//3 = and the conclusion is the same. 

We can try to take the limit of the fundamental solution that exists for m > m c and 
prove that the solution concentrates around the origin as m — > m c . This is clear in 
the standard FDE since the solutions are explicit, see the very peculiar case N = 2, 
m — > in Lemma 8.3. 



An analysis of the behaviour of the potentials is also illuminating. We leave it to the 
reader. 
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12 Precise decay rates of the Barenblatt profiles 



We resume the study of the decay rates of the Barenblatt solutions started at the end 
of Section [8] 

• Case m > 1. We may get the exact rate as follows: since F m < Cr~ mN we 
obtain by use of the convolution formula that L s iF m ~ r -N+2s ^ an( j means ^hat 
f °° rF(r)dr behaves like r~ N+2 ~ 2s as r — > oo. But this also implies that the "mass" in 
an annulus given by the integral J^ r rF(r)dr behaves like r~ N+2 ~ 2s (with a different 
constant). Using the monotonicity of F we arrive at the consequence that F has the 
decay rate 



fl2.1) F(r)~Cr 



-N-2s 



valid for all m > 1, which equals the decay rate for m — 1. This decay for s < 1 is in 
stark contrast with the case s — 1, where the profiles of the Barenblatt solutions are 
compactly supported in the same m-interval, m > 1. 

• Case m c < m < 1. A power-like bound from below is obtained as follows: We may 
start from the homogeneity estimate p2] that says that (1 — m)tut < u. In terms of 
the self-similar profile, this just means that —(1 — m)f3(N F + rF'(r)) < F, hence 

-rF'(r) AT 1 2s 

<N + 



F(r) (l — m)/3 1 — m 

Integration of this inequality gives the following lower bound, valid for all r > 1, all 
s G (0, 1) and all m > m c (s, N): 

(12.2) F(r)>Cr- 2s/(1 - m) . 

Moreover, the function J(r) := F(r) r - 2s K l - m ) i s monotone non-decreasing with r, so 



that it has a limit as r — > oo. The lower bound (12.2) is a good starting point since it 
is the exact decay rate for the standard diffusion case s = 1 when m c < m < 1 with 
compactly supported initial data, 



We will obtain a similar upper bound in the range m c < m < mi as a consequence of 



the VSS construction in Section 13, where the asymptotic constant c^ = lim^oo J(r) 



is calculated. See more details in Section 13. Note in passing that 2s/ (1 — m) > N 



precisely for m > m c (s,N). Summing up, the decay rate (12.2) is optimal for m < 
mi = N/(N + 2), and a very precise rate is obtained. 
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• However, the rate (12.2) is very far from a realistic estimate for m close to 1, since 
2s/(l — m) — > oo, which is not, loosely speaking, an admissible decay in the situation 
of fractional Laplacian diffusion. 

We look for a bound from above as follows. Since F(r) < C r~ N we have F m < 
Cr -Nm and the fact that Nm < N imp i ies that L s /F m < Cr~ Nm+2s \ hence 

/oo 
rF{r)dr < C r ~ Nm+2s ' . 

Using the monotonicity of F this means that F(r) < C r - mN - 2s . But Nm + 2s > N 
for m > m c , so that we have a gain of upper bound exponent from 70 = iV to 
7x = mN + 2s. Proceeding iteratively as long as the whole argument we have used 
is justifed, we get a sequence of increasing exponents 7^ given by j k+1 = j k m + 2s, 
that converge to the fixed point of the iteration formula, i. e., 2s/ (1 — m). This is all 
right for m < ni\. However, for 1 > m > ni\ we have 2sm/ (1 — m) > iV so that after 
a number of steps F(r) m < r~ mik with 7^771 > N. In that case, the next line of the 
previous argument changes to give L s /F m ~ r -N+2s ^ an ^ following the same line of 
argument F(r) ~ t--^- 2 ^ for all large r,which is the desired conclusion. 

• In the limit case m = N/(N + 2s), the iteration from above is not inrerrupted 
and we get F(r) < C r~ N ~ 2s+£ for every e > 0. Regarding the iteration from below, 
the starting rate is 7 = iV + 2s, so that the exponent of our lower bound for F m is 
mN/ (iV+2s) = N . In this case the convolution formula for L s > produces a logarithmic 
correction for L T iF m that we have to take into account. The estimate is then 

(12.3) F{r) > cr~ N ~ 2s log(r), for all large r; m = m\. 
We may sum up the results as follows. 

Theorem 12.1 For every m > mi = N/(N + 2s) we have the asymptotic estimate 

(12.4) lim F M {r) r N+2s = d M a , 

r— >oo 

where M = J F(x) dx, C\ = Ci(m, N, s) > 0, and a = (m — mi)(n + 2s)/3. On the 
other hand, for m c < m < m\, there is a constant Cooim, N, s) such that 

(12.5) lim F M {r) r 2s ^ 1 -'^ = C^. 

The case m = mi is borderline and has a logarithmic correction. 
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The analysis done before assumed M — 1, for M ^ 1 just use scaling. For m G 
(m c ,mi) the fact that does not depend on M is commented and explained in 
whole detail in Section 13 where the last elements needed to complete this proof are 
given. 




Figure 2: Plot of decay rates of the self-similar profiles. Here, N = 3 and s = 0.8. 



13 Very singular solutions in the fast diffusion range 



13.1. Existence. A special solution that plays an important role in the theory of 
the standard Fast Diffusion Equation is the so-called Very Singular Solution (VSS), 
cf . [231 EH] . More precisely, foer every 1 > m > m c there is a global in time function 
with the form of separation of variables that solves the FDE away from x = and 
has a standing singularity at x = for all times. It also starts with V(x, 0) = for 
x ^ 0. 

Using the same idea, we want to find a similar solution for the fractional diffusion 
equation for < s < 1 in the good fast diffusion range m c < m < 1. We look for 
solutions of the form 

U(x,t) =T{t)X{x). 
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Substitution into the equation leads to the value T(t) = t 1 ^ 1 m ) for the time factor, 
while Y = X m has to satisfy 



(13.1) 



-A) S Y 



l — m 



Y p = 0, p= 1/m. 



Now we try as solution the function Y(x) = C m \x\ a . Under suitable conditions on 
a, to be discussed below, we will have 



(13.2) 



_ A yy = C m k\x\- {a+2s) . 



for some negative constant k = k(a,N,s). It follows that Equation (13.1) for Y can 
be satisfied only if a takes the value a(m,s) = 2s/ (p — 1) = 2sm/{l — m). The 
constant c > is then determined by —C' m k(a) = C/(l — m), so that 



(13.3) 



C 



l—m 



m)(—k(a)). 



Let us examine the conditions on a. 



(i) First, we need the condition a < N to make Y integrable and formula (13.2) hold 



which means m < mi := N/(N + 2s). Without this restriction we enter into a theory 
of Laplacians of non-integrable functions, which is acceptable anyway for s = 1, but 
it does not work in the fractional case s < 1, as we will see below. 



(ii) We need to calculate k to know when it is negative. This is done in (15.3). Note 
that in the interval m c < m < mi we have iV — 2s < a < N so that N < a + 2s < 
N + 2s. We check there that k(a) > for all a > N - 2s. 

Finally, let us check the properties of the constructed function. We remark that U 
is not a weak solution in the standard sense since U(-,t) is not even integrable in 
space around the origin. However, U m is locally integrable, and we can propose a 
generalized definition of weak solution that restricts test functions to be supported 
away from the origin. We denote this situation with the name generalized weak 
solution with an isolated singularity. 



Theorem 13.1 There exists a VSS for Equation (15.1) in the range m c < m < mi 
and it has the explicit form 



(13.4) 



U(x,t) = Ct 1/{1 - m) |x|- 2s/(1 ~ m) 



where C = C(m,s,N) is given by the explicit formula (13.3). U is a weak solution 
in a generalized sense, more precisely, a weak solution with an isolated singularity. 
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In order to compute the exact value of C we must use (13.3) together with (15.3). 
For further reference we callthe constant Cyss- 



1131 2. VSS as limit. We now establish the relation of the VSS with the Barenblatt 
solutions. 



Theorem 13.2 The VSS (13.4) is the limit of the Barenblatt solutions u* M (x,t) as 
the mass M goes to infinity. 



Proof. We need to prove that this VSS is the limit of the Barenblatt solutions as the 
mass goes to infinity. By the known comparison properties the sequence of Barenblatt 
solutions u* M (x,t) is monotone increasing with M > 0. 

Next, we check that they are all bounded above by the VSS. This is done by compar- 
ison. But instead of comparing u* M with U, we compare u(x,t) and U, where u is a 
solution with mass M and smooth data Uq supported in a small ball, so that u(x, 0) < 
U (x, 1) for all x e K . The comparison result holds and proves that u(x, t) < U (x, t) 
for all t > 0. Now we do the rescaling that is used in the proof of asymptotic behaviour 
and observe that U is invariant. Therefore, T\u(x,t) < U(x,t). In the limit t — > oo 
we use the result on asymptotic convergence to get u* M (x, 1) < U(x, 1). The result is 
obviously true for all other times with a minimal modification of the argument. 

This upper bound allows us to pass to the monotone limit in the family {u* m (x,i) : 
M > 0} and obtain a function U*(x,t) < U(x,t). This upper bound implies that U* 
is a weak solution of the Fractional FDE, and it is locally bounded away from the 
space origin x = for all times t > 0. We conclude that U* is another possible very 
singular solution. 

It also has the same invariance under (T^u)(x,t) = yU 2s ^ 1_m )n(/is, t) which means 
that the space form must be c(i:)|a:| 2s//( ' 1 ~ m ^. But this implies the self-similar form we 
have chosen and this implies that U* = U 

|13[ 3. Precise asymptotic behaviour. We claim that the asymptotic constant 

(13.5) Coo (F M ) = lim F M (r)r 2s ^-^ 

for the Barenblatt solutions coincides with the constant C for VSS, let us call it Cyss- 

Indeed, by the ordering of the solutions we must have c^Fm) < Cyss- On the 
other hand, the scaling of the Barenblatt solutions immediately implies that c^Fm) 
is independent of M. Therefore, for every M > we have 

(13.6) F M {r)< Coo r- 2s ^- m \ 
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If we had Coo(F) < Cyss, then the limit of the Barenblatt solutions as M — > oo would 
be a very singular solution of the form (13.4) with a constant less than Cyss, and 
this is not possible, as we have already seen. 



13,4. Infinite limit in other ranges. There are no fundamental solutions for 



m < m c so that the ranges to discuss are m > 1 and m\ < m < 1. In the linear case 
m =! the scaling of the Barenblatt solutions is 

(13.7) u* M (x,t) = Mu(x,t) , 

so it is clear that limjv/^oo u* M (x, t) = oo everywhere in Q = R N x (0, oo). A similar 
result holds for m > 1, in view of the scaling 

(13.8) u* M (x, t) = ^ 2s u{^- m x, t), M — ^{m-l^s _ 

Since u is monotone in r = \x\ and m > 1, we have u* M (x,t) > fi 2s u(x,t) and the 
result follows letting M — > oo. 

In the remaining fast diffusion range, m x < m < 1, the final answer is the same, 
an infinite limit, but the argument is not so simple. It is based on the minimum 
decay rate of the Barenblatt profiles, like 0(r~( N+2s ^) (with logarithmic correction if 
m = mi). In particular, such behaviour implies that the limit 

lim F 1 (r)t 2s/(1 - m) = +oo. 
r— >oo 

The conclusion follows by scaling: For all large M > and putting M = ^ N ( m - 1 )+ 2s 
we get for fixed r > as M — >■ oo: 

lim F M (r) = lim fj, 2s F^^r) = r ~ 2s/{1 - m) lim z^F^z) = oo. 



Remarks. (1) This last divergence result means that no VSS in the sense of the 
beginning of this section can exist in this range. The result represents a marked 
difference with the case s = 1, where a VSS solution exists for all m e (m c , 1), and 
the Barenblatt solutions are uniformly bounded by it. 

(2) The construction of a separated-variables solution with a fixed isolated singularity 
can de performed for m < m c but then we would get finite-time extinction profiles, 
in the line of |55j, Chapter 5 . 
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14 Construction of eternal solutions for m = m c 

We have already seen that the Barenblatt solutions in the sense of fundamental solu- 
tions with self-similar form, do not exist for m < m c := (N — 2s) /N. One wonders if 
there is a special way to pass to the limit m — >■ m c to obtain a self-similar solution, 
even if it cannot be a fundamental solution. 



Theorem 14.1 There exists a family of eternal solutions of equation (1.1) with ex- 
ponent m = m c and < s < 1 . They have the form 

(14.1) U(x,t) = e- Nct F{\x\e- ct ) 

where F(r) is a continuous, radial function with F(0) = a > and F(r) r N — > b > 
as r — > oo. 

Proof. 1. The process is known in the standard case s — 1, but it will be convenient 
to explain it in some detail in this context, since it involves a subtle manipulation of 
the equations. We consider a different scaling of the flow that is good when passing 
to the limit m — > m c = (N — 2)/N with m > m c . It says 

u(x, t) = (1 + (ct/f3))- n ^v(y, r), y = x(l + (ct/f3))^, r = (3 log(l + c(t/f3)) . 

It applies for t > — c/3, where /3 is the self-similar exponent and c > is a free 
constant. The original equation]^] u t = A(u m /m) becomes 

(14.2) cv T = A(v m /m) + cV • (yv) = V(vV(-p+ -y 2 )), p=^—v m -\ 

2 1 — m 

Notice that no /3 appears, which is important since (3 — > oo as m — > m c . Equation 



(14.2) has the stationary solution (profile) 



with b > arbitrary. This holds for 1 > m > m c . In the limit m — > m c we have: 
1/(1 -m) -)• N/2, (3 -> oo, r(t) -> t, (1 + (ct//3))^ -> e" rf , (1 + (ct//3))~ N ^ -> e" cAr *, 
so that the profile in the limit is 

^ c (y) = (&+^ 2 r Ar/2 , 



3 Note the slight change in the form of the equation due to the constant m, that is custom- 
ary in studies of FDE. It is inessential here but we have inserted it for agreement with published 
information. 
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and the solution becomes (with a = c/N) 



(14.3) U mc (x,t) = e- cNt (b + ax 2 e- 2ct )- N ' 2 = {ax 2 + be 2Nat )- N ' 2 , 

which is the (two-parameter family of) eternal solutions obtained in [55], formula 
(5.36). We may still enlarge the family by displacing the origin of space coordinates. 
The scaling group Tu(x,t) = X N fx~ N ^ 2 u(\x, /it), A, fj, > 0, acts on these solutions to 
change the parameters a and b in the form: T(a, b) = (a/z, bfi/\ 2 ), so the whole family 
of solutions comes from a single one. 

Moreover, when b — > we get the singular stationary case [/^(x, t) = A\x\ , 
A = a~ N l 2 > 0, which is not a solution at x = since A£7™ has a delta function 
there. This situation reminds of the increasing limit of Barenblatt solutions to the 
VSS for m > m c , in that case we had used the mass as paramter, here the parameter 
is rather F mc (0) = b~ N l 2 and a > is kept fixed. Note also that does not depend 
on time, and that the constant A is arbitrary. 

2. We now apply these ideas to the limit m — > m c (s) = (N — 2s) /N in the fractional 
diffusion case < s < 1. This case is more difficult since the Barenblatt solutions are 
not explicit. We have to avoid the vanishing of the space asymptotic constant in the 
limit by careful time scaling of the above type that leads to the correct renormalized 
equation. 

Indeed, the equations for m> m c have VSS solutions in the standard scaling of the 
form 

V(x) = ct l/{l ~ m) x~ 2s/{1 - m \ 
where Cyss is given by C x ~ m = — (1 — m)k(a) and 



k{a) = 2 



2s T((N-a)/2)T((a + 2s)/2) 
r((JV-a-2s)/2)r(a/2) 



with a = 2ms /(l — m). When m — > m c we have a ~ iV — 2s, (N — a)/2 ~ s, 
(a + 2s) j2 ~ N/2, a/2 ~ (JV/2) - s, and finally 

N - a - 2s N -2s sm N(l - m) - 2s 1 



m 2(1 -m) 2(l-m)/3 



Therefore, C l ~ m ~ Comparing the profile equations (8.1) and the one above, 

we see that we replace (3 by c so that the new constant C is given by (7 1_m = k\/c 
which can be made 1 by suitable choice of c. 
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3. Construction. We now consider the family of self-similar solutions for m ~ m c 
with F m (0) = 1 (which means that you need to choose the mass of the solution 
appropriately). We also have 

F m (r) < C m r- 2s ^- m \ lim F m {r)r 2s '^ = C m 1. 

Using the monotonicity of the profiles and C a regularity of the family F m we may 
pass to the limit F m (r) -> F(r) so that F(0) = 1, F(r) < r~ N , with F G C a , and 
F(r) r N increasing. This is the eternal solution, after some extra work that is rather 
easy. 

You may then construct a two parameter family of radial solutions, and then add 
displacement in space. The two parameters are F(0) controlled by constant mass 
scaling, and the constant at infinity that is controlled by the other part of the scaling 
group. □ 



15 Aleksandrov's symmetry principle 

The Aleksandrov-Serrin reflection method is a well-established tool to prove mono- 
tonicity of solutions of wide classes of (possibly nonlinear) elliptic and parabolic equa- 
tions, cf . [31 EI] ■ It has been quite useful in particular in the case of the PME, as 
documented in [56J. Here we will establish a version of the principle valid in the 
presence of fractional operators, and then derive useful monotonicity results. 



1. Elliptic setting. We consider the problem 
(15.1) Lu m + u = f 

posed in W N with L = (— A) s , m > 0, and / G L 1 (IR Ar ). We take a hyperplane H that 
divides 1R into two half-spaces Qi and Q2 and consider the symmetry with respect 
to II that maps Qi onto f^. Then we have. 



Theorem 15.1 Let u be the unique solution of (15.1) with data f G L 1 (IR Ar ) ; / > 0. 
Under the assumption that 

(15.2) f(x)>f(U(x)) in ^ 
we have 

(15.3) u(x) > u{U{x)) in fi x . 
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Proof, (i) Due to the translation and rotation invariance of the equation we may 
assume that H = {x G R N : Xi = 0}, and fli = {x G R N : x x > 0}. We write 
x' = n(x) where II is the symmetry with respect to x 1 = 0. We take a function / 
such that f(x') < f(x) when Xi > 0. We assume at this stage that / is bounded, 
continuous and integrable. 

We solve the elliptic problem with data / to get a solution u and write u(x) = u(x'), 
that solves the problem with data f(x) = f(x') in the whole space R N . 

(ii) We now assune for a moment that u < u in Qi and consider a point xq G Oi 
where u touches u from below. If such a point exists, then we have u(x') < u(x) for 
all x G Vt\ and we get that 



Lu m (x ) = C J 



u m {x Q )-u m {y) 



Fo - y 



N+2s 



dy 



and 

an*) - c i ^ - gi (y> dv= c( "7") > 

Therefore, when calculating Lw m (a;o) — Lu m (xo) the terms in -u m (a;o) = ^(xq) cancel 
out and we get 

is-teO - - c / ""("'-y «, + c [ rf9 . 

Aii Fo-j/r +2a Jn 2 Fo-yr +2s 

Interchange of y and ?/ in the last integral gives 

r-m, \ t m t \ f u m (y)-u m (y') , f u m (y')-u m (y) , 

Lu m (x )-Lu m (x )= i m£r d y+ i 7i7v+2s y ' 

hence, 

L(Q» - »-)(«,) = f^y) - »- M " R^Til) * > 0- 

where we have used that for x , y G f2i we have |x — y| < \x — y'\. On the other hand, 
at this point we have f(xo) > f{x' ). We get a contradiction with the equation unless 
both solutions coincide everywhere. The solution must then equal its symmetric 
reflectio, which is impossible is / is not symmetric. 

(iii) In order to apply this argument we need to show that the situation where the 
solution and its reflection are ordered happens, which is not clear. The general 
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argument that fits the actual situation is done by arguing on an approximation of 
the solution, as typical in the proofs of maximum principles for standard elliptic 
equations. We argue as follows: 

- We take / in, say, L°°{R N ) nL 1 (M JV ) and with compact support. We use the theory 
of existence for data / = / + C, which is not difficult since it is like the one done 
in the papers [12J H6] after a vertical displacement of the solution. We thus get for 
bounded data fc = f + C, C > 0, a solution uc > C and uc — C G L 1 (WL N ). 

- Under the stated hypotheses the existence and regularity theory says that uc G 
L°°(M. N ) n L 1 (IR Ar ), is continuous and that uc > uo = u. The regularity is obtained 
by a Morrey embedding theorem. 

- A known strong maximum principle for fractional Laplacian operators implies that 
these functions do not touch in M. N , i. e, uc > u . The argument is not very different 
from Step (ii) above, but simpler since there is no reflection. 

- Let us now check the symmetrization argument on this family. We will write for 
more clarity uc = u(f + C), u = u(f). For large C > we have u(f) < u(f + C) 
and u(f) < u(f + C) in Qi («(/) indicates the symmetrical image of u(f) as before). 
We use the argument displayed before to conclude that there is no possible point of 
contact for u(f) and u(f + C) inside Q\ and there was none at the boundary H where 
<f) = u(f). 

- A this point we may lower the C until we get a point of contact in one of the two 
comparisons. If C > there is no contact at X\ — and none as \x\ — > oo and we get 
a contradiction in the two comparisons. 

Therefore, the infimum of the C is zero. Then, there is contact as \x\ — > oo and at 
the boundary Xi = 0. But the plain comparison u(f) < u(f). The proof is finished 
under the stated assumptions on /. 

(iv) For general / G L 1 (IR Ar ) we prove the result by approximation with functions f n 
as above and then use the /^-continuity of the map / h- > u(f), jlS]- □ 

II. Parabolic setting. We consider the problem 
(15.4) u t + Lu m = 0, 

posed for x G WL N and t > with initial data u(x, 0) = u (x) G L 1 (IR Ar ), u > 0. 
We use the same notations and conventions for the hyperplane H, half-space Qi and 
symmetry II. Then we have 
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Theorem 15.2 Let u the unique solution of (15.4) with initial data uq. Under the 
assumption that 



(15.5) Uq(x) > uo(H(x)) in Qi 
we have for all t > 

(15.6) u{x,t) >u(U{x),t) for x e tti. 

Proof. This is an easy consequence of the fact that the nonlinear operator S : u t— > 
Lu m can be defined as m-accretive in / 1 (IR iV ), and its resolvent (/ + AS 1 ) -1 also and 
ordered contraction, cf. [16]. We can then solve the parabolic problem by implicit 
discretization in time, i.e., for times to — 0,£i = h,t 2 = 2h,--- , we define the 
approximations 

^(u(t n -u(t n _ 1 )) + L(u(t n ) m ) = 0. 

We apply the conclusion of the elliptic theorem iteratively at every step to get 
u(x,t n ) > u(Tl(x),t n ) for every x G n > 1. Finally, we pass to the limit h — > 
in the discretization step (using the Crandall-Liggett Theorem as in [46J) to get 
the solution u(x,t) with the desired properties. See more details of this method in 

[251 EH EE]. □ 

III. Application. Comparison after reflection around suitable hyperplanes allows 
now to conclude that the solution with initial data supported in a small ball B £ will 
be decreasing in the radial exterior direction in the annulus {x : \x\ > 2e}. We also 
conclude that it is decreasing along a cone of directions centered along the radius 
with aperture that goes to the flat cone when the annulus is taken {x : \x\ > R} with 
R/e oo. This implies that the solution with initial date a Dirac delta must be 
radially symmetric. All these techniques are well-known, as explained in Chapter 9 
of [56], where references to further literature are given. 

Appendix 1. Fractional Laplacians and potentials 

According to Stein, [52], chapter V, the definition of (— A)^/ 2 is done by means of 
Fourier series 

(15.1) ((-A)^f)(x) = (2n\x\ff(x) 
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and can be used for positive and negative values of f3. For f3 = —a negative, with 
< a < N, we have the equivalence with the Riesz potentials 



(15-2) (-A)-/ 2 / = I a {f) := f | fiV L a dy 

7(a) J r n \x-y\ N a 

(acting on functions of the class S for instance) with precise constant 

7 (a) = 7r N/2 2 a T(a/2)/T((N - a)/2). 

Note that 7(a) — > 00 as a — > N, but j(a)/(N — a) converges to a nonzero constant, 
7r N / 2 2 N ~ 1 T(N/2). Note that in this notation a stands for 2s in our previous notation. 

• The Fourier Transform of the function f(x) = \x\~ N+a is /(£) = 7(a) (2^)"" |^| 
The Fourier Transform of the function f(x) = \x\- a is /(£) = 7(A^-a)(27r) a - Ar |^| a - iV . 
We conclude that its s-Laplacian, (— A) s \x\~ a , is the power k(a) |:r|~ a ~ 2s , with a 
constant factor that equals 

7(iV-a) 7(iV-a) 



fc(a) 



7(A^-a-2s) 7(A^-a-2s) 
so that 



2 ,r((iV-a)/2)r((a + 2 S )/2) 



(15 - 3) Ka)=2 T((N-a-2s)/2)T(a/2) 

When s = 1 we get 

k(a) =a(N -a-2) 

which is positive whenever a G (0, N — 2). For < s < 1 it is not so clear. In any case, 
in the range of interest, < iV — 2s < a < n, all the factors minus one are positive and 
the remaining one, T((N — a — 2s)/ 2) corresponds to an argument ((N — 2s) — a)/2 
which is negative but larger than — s > —1, hence T((N — a — 2s)/2) < 0. 

Appendix 2. Bessel kernels 

They are usually introduced via Fourier transform, 

G a (0 = (l + |£lT a/2 , £eR N ,a>0. 
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There is an expression for the kernel of the form 

roc j. 

G a (x) = c{a) / e^l 2 /' &~ N ^ 2 -, xeR N , 
Jo t 

see Section V.3.1, so that G a is a non-negative, radially decreasing function. 
Moreover, G a is integrable and c(a) is a positive constant chosen so that = 1. 

The Bessel potential of order a > of the density p is defined by 

B a (x) = / G a (x -y)p{y)dy. 

We have the following estimates for < a < N: 

< G Q (x) < C \x\- {N ~ a \ if < |x| < 1, 

< G a (x) < e" |x|/2 if \x\ > 1 , 

where C = C(a;N). More precisely, the kernel can be represented by means of the 
McDonald function: 

G a {x) = c(a,N)\x\^~ N ^ 2 K {N ^ a)/2 (\x\). 

The Macdonald function with index is, K v , is 6 R, called also the modified Bessel 
function of the second kind, is given by the following formula: 

/■oo 

K v {r) = 2- l - u r u / e-^^r^dt, r > 0. 
Jo 

The asymptotic behaviour of K v is as follows: 

K A r ) ~ ^Y^( r / 2 ) _iy as r -> 0+, i/ > 
while Ko(r) ~ log(l/r) as r — )■ 0+. At infinity we have 




as r — t- oo, 



where the notation g(r) ~ f(r) means that the ratio of g and / tends to 1. For v < 
we have K_ v (r) = K u (r), which determines the asymptotic behaviour for negative 
indices. For properties of K v see [27]. 
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Comments and extensions 



• The linear case was known, and our results recover the information with a completely 
different machinery. 

• One wonders if there are any explicit or semi-explicit formulas for the family of 
Barenblatt solutions of this paper. 

• The limits s — > and s — y 1 are worth examining. 

• Fundamental solutions can be constructed for the equation posed in a bounded 
domain. They do not play however such an important role in the theory. For instance, 
they shed no light on the asymptotic behaviour, even in the simplest case of zero 
Dirichlet boundary conditions. 

• The positivity of the fundamental solutions in all M , and more precisely, the 
power-like maximal decay rate as \x\ — > oo, are properties shared by all nonnegative 
solutions of the FPME and not only the solutions that take a special self-similar 
form. Quantitative versions of the lower bounds on the positivity of solutions are 
established in a paper with Bonforte [18]. Higher regularity of positive solutions of 



equation (1.1) is currently being investigated. 

• There is the open problem of proving the uniqueness of solutions with measure 
data, not just Dirac deltas. 

• In the case of the standard PME, the range m < m c is now understood: no 
fundamental solutions exist, but the family of Barenblatt solutions for m > m c can 
be continued algebraically and they exhibit interesting attraction properties, [T5l [T7] . 
Not much is known to our knowledge for the corresponding fractional equations. 

• New applications of these models in the applied sciences would be welcome as a 
motivation to develop further aspects of this theory. 
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